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ABSTRACT 

A full set of calculations is presented for inner-shell n = 2 to 3 photoexcita- 
tion of the 16 iron charge states: Fe 0+ (Fe I) through Fe 15+ (Fe XVI). The blend 
of the numerous absorption lines arising from these excitations (mainly 2p - 3d) 
forms an unresolved transition array (UTA), which has been recently identified as 
a prominent feature between 16 - 17 A in the soft X-ray spectra of active galactic 
nuclei (AGN). Despite the blending within charge-states, the ample separation 
between the individual-ion features enables precise diagnostics of the ionization 
t> ; range in the absorbing medium. Column density and turbulent velocity diagnos- 

tics are also possible, albeit to a lesser accuracy. An abbreviated set of atomic 
parameters useful for modeling the Fe 2p - 3d UTA is given. It is shown that the 
O ' effects of accompanying photoexcitation to higher levels (n > 4), as well as the 

associated photoionization edges, may also be relevant to AGN spectra. 

Plh| Subject headings: atomic data — atomic processes — line: formation — galaxies: 

active — galaxies: nuclei — X-rays: galaxies 

ci 

1. INTRODUCTION 

Active galactic nuclei (AGN) viewed directly towards the supermassive black hole (e.g. 
Seyfert 1 galaxies) can produce a rich absorption spectrum. High-flux X-ray and UV contin- 
uum emission from the central source ionizes the ambient gas to high charge states. Conse- 
quently, the ionized gas along the line of sight imposes absorption features on the observed 
continuum. Away from the central source, the ionizing flux gradually diminishes, limiting the 
maximum charge states that can be formed, especially for high-Z elements. Notwithstanding, 
the photon energy distribution, generally a power law, allows for a substantial fraction of 
energetic photons. Therefore, the intermediately ionized ions are exposed to X-ray photons 
that can be appreciably more energetic than the typical energies for ionization or excitation 
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from the valence shells of these ions. When the column density through this gas is sufficient, 
absorption by means of inner-shell photoexcitation processes plays an important role in the 
formation of discrete features in the spectrum. 

In particular, the n = 2 to 3 absorption lines (n being the principal quantum number 
of the active electron) of the cosmically abundant element iron (Fe) fall in the soft X- 
ray band that can be readily measured by the grating spectrometers on board Chandra 
and XMM-Newton. The n = 2 to 3 valence transitions are well known in the soft X- 
ray band for the (L-shell) ions Fe 16+ to Fe 23+ . The analogous, but inner-shell, excitations 
in Fe 0+ to Fe 15+ (referred to hereafter in shorthand as Fe M-shell ions), have not been 
given the same attention. In the laboratory, Chenais-Popovics et al. (2001) studied the 
opacity of a thin, laser heated Fe foil illuminated by a strong X-ray backlighter. They 
detected unresolved inner-shell n = 2 to 3 photoabsorption by Fe 4+ to Fe 9+ around 17 A and 
modeled it assuming local thermodynamic equilibrium. Recently, we have reported the first 
astrophysical detection of a similar feature in the spectrum of the luminous quasar IRAS 
13349+2438 (Sako et al. 2001) obtained with the Reflection Grating Spectrometer (RGS) 
on board XMM-Newton. In that spectrum, we ascribe the broad absorption trough between 
16 - 17 A to an unresolved transition array (UTA) of inner-shell n = 2 to 3 (mainly 2p - 3d) 
lines pertaining to various Fe ions. As shown in Sako et al. (2001), the exact shape of this 
absorption feature is well fitted with a superposition of lines from Fe 6+ to Fe 11+ with ionic 
column densities of 0.9 - 6.4 x 10 16 cm" 2 , which correspond roughly to column densities of 
hydrogen atoms of a few times 10 21 cm -2 . Recently, the soft X-ray spectrum of the Seyfert 1 
NGC 3783 obtained with the RGS has revealed a similar, but even more prominent Fe 2p - 
3d UTA (Blustin et al. 2001). This feature, unaccounted for in the analysis of the Chandra 
spectrum of NGC 3783 (Kaspi et al. 2000; Kaspi et al. 2001), could explain the discrepancies 
between the data and the model around 17 A in that work. A much weaker Fe 2p - 3d UTA 
was tentatively identified in the Seyfert 1 Markarian 509 by Pounds et al. (2001). Fig. 1 
shows the Fe 2p - 3d UTA region of the RGS spectrum for both IRAS 13349+2438 and 
NGC 3783. Despite the apparent difference in absolute flux, the shape and structure of 
the UTA are somewhat similar, covering approximately the same range of ionization. The 
transmission in NGC 3783 is evidently lower. 

For both IRAS 13349+2438 and NGC 3783, as well as for other AGN, the Fe 2p - 
3d UTA can be used to diagnose the ionization state and column density for the relatively 
low-ionization components of the soft X-ray absorbing gas. Since the Fe M-shell ions also 
absorb in the UV and EUV bands, they can be used to relate the UV absorber with the 
X-ray absorber, which have curiously shown different column densities for many sources. 
Additionally, the same charge states that produce the Fe 2p - 3d UTA absorption also 
produce strong emission lines in the various optical bands. If both features are observed 
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simultaneously, the Fe 2p - 3d UTA can potentially provide an independent probe of the 
level of ionization in the optically emitting gas. Furthermore, if the low-ionization Fe M- 
shell gas is associated with the ionized skin of the putative AGN torus, as postulated by Sako 
et al. (2001), it can provide a constraint on the size of the torus and the inclination angle 
of the (host) galactic plane with respect to our line of sight. 

Photoabsorption from the L-shell (n= 2) of neutral (or nearly neutral) Fe falls between 
17 and 18 A. This was most recently detected by Paerels et al. (2001) in the low-mass 
X-ray binary X0614+091. Paerels et al. (2001) mention only the L-shell photoionization 
edges, probably because of the lack of atomic data for the associated (photoexcitation) 
lines. Owing to their much higher cross sections, lines should always appear when edges are 
detected. Thus, in addition to AGN observations, the present calculations are also relevant 
to spectra produced by absorption of local X-ray sources, such as X06 14+091, by interstellar 
Fe in our galaxy. This paper provides an account of the atomic processes and the full details 
of the spectral modeling issues related to the inner-shell n = 2 to 3 absorption lines of Fe 
M-shell ions. A compact set of wavelengths and oscillator strengths is made available for 
fitting this feature to observed spectra and for diagnostics. 



2. ATOMIC CALCULATIONS 

In order to calculate the wavelengths and oscillator strengths for the radiative tran- 
sitions, we use the multi-configuration, relativistic HULLAC (Hebrew University Lawrence 
Livermore Atomic Code) computer package (Bar-Shalom, Klapisch, & Oreg 2001). The 
intermediate-coupling level energies are calculated using the relativistic version (Luc-Koenig 
1972; Klapisch et al. 1977) of the Parametric Potential method by Klapisch (1971). The 
relativistic approach employed in HULLAC, which uses the jj coupling scheme for its initial 
vector basis, is perfectly adequate for the inner-shell excitations studied in this work. The 
high efficiency of HULLAC is also crucial, owing to the numerous transitions and the consid- 
erable complexity of the electronic structure in the ions considered. Calculations are carried 
out for neutral iron (Fe 0+ ) through Na-like Fe 15+ . The lower Fe charge states could be less 
relevant to the AGN vicinity, but more to the narrow line regions or to X-ray illuminated 
interstellar gas. 

The atomic structure calculations include the ground configurations 2s 2 2p 6 3P (x = 1 
to 14) for Fe 15+ through Fe 2+ , and 2s 2 2p 6 31 14 4s :c (x = 1, 2) for Fe 1+ and Fe 0+ . For the 
ground configurations of Fe 1+ and Fe 0+ , the [3d + 4s] configuration mixings are included, 
i.e., 2s 2 2p 5 3s 2 3p 6 [3d 6 4s x + Sd^s*" 1 ] (x = 1, 2). All of the n = 2 to 3 photoexcitation 
transitions from the ground level are calculated for each ion, among which the 2p-3d and 
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2s-3p are the strongest. For charge states Fe 8+ and lower, the 3p subshell is full and, thus, 
the 2s-3p excitation channel is not possible. For Fe 15+ , the 2p-3s excitation channel is open 
and is also calculated. For the 2p-excited Fe 10+ - Fe 15+ ions, the important [3s3d + 3p 2 ] 
configuration mixings are possible and are included, i.e., 2s 2 2p 5 [3s 2 3p :c 3d + 3s3p :c+2 ] (x — 
to 4), as well as 2s 2 2p 5 [3s3d + 3p 2 ] for Fe 15+ . The significance of other mixing effects, such 
as [3p 2 + 3d 2 ], which were found important for excitations out of the 3p subshell (Quinet 
& Hansen 1995), were also examined. We find these mixings to reduce the total n = 2 
to 3 photoexcitation effect by about 1% and to alter the level energies by much less than 
that, while increasing the number of absorption lines by a factor of ten and more. Therefore, 
for the sake of simplifying the computations, we have excluded these mixings. Electric and 
magnetic dipole and quadrupole transitions were computed. The contribution, however, of 
excitations other than electric-dipole is less than 0.1 %. 

The inner-shell excited levels studied here tend to strongly autoionize. This has several 
effects on the observed spectrum; First, through the low fluorescence yield. Since the inner- 
shell n = 2 to 3 photoexcitation processes are by and large followed by autoionization rather 
than X-ray re-emission, they leave an observable imprint only in absorption, but not in 
emission. This absence of fluorescence following L-shell photoexcitation is evident in the 
soft X-ray spectra of Seyfert 1 galaxies, but even more so in the Seyfert 2 spectra, where 
emission is not confused with absorption. The low fluorescence yield further implies that 
dielectronic recombination via these channels is inefficient. Another consequence of the 
high autoionization rates is the broad natural widths of the absorption lines associated 
with inner-shell photoexcitation. At high column densities, the natural width of the line 
can become important in determining the total flux absorbed. For many of the transitions 
considered in this work, the total autoionization rate is an order of magnitude higher than 
the radiative decay rate. Thus, in order to obtain correct line profiles and to account for the 
total photoabsorption effect, we have calculated the total depletion rate incorporating both 
the autoionization and radiative decay rates to obtain the natural line width. The natural 
widths obtained here could be, in some cases, a slight underestimation of the actual width 
because including all of the autoionization transitions is computationally not permissible. 
The ionized configurations that were included in the autoionization calculations were the 
most important ones that can be reached by a single-electron transition from the ground 
level. 

The results of the above-described calculations are shown in Figs. 2 and 3, where the 
absorption spectra produced by the inner-shell n = 2 to 3 excitations are presented for Fe 0+ 
through Fe 15+ . The complex atomic structure and the vast amount of detail in the spectra 
can be clearly seen in the figures. For display purposes, a low ionic column density of Ni on = 
2 xlO 15 cm -2 and a temperature of kT ion = 10 eV are chosen, with no turbulent velocity. 
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These parameter values help minimize the effects of saturation and blending of the individual 
spectral lines in the plot. The dependence of the observed spectrum on these parameters is 
discussed below in more detail. 

The largest group of lines for each ion represents the 2p - 3d transition array, and can 
be seen in the figures to shift from ~ 17.5 A for Fe 0+ to ~ 15.2 A for Fe 15+ . In fact, the 
predominant transitions are 2p 3 / 2 - 3d 5 / 2 and 2p^ 2 - 3d 3 / 2 • The 2p 3 / 2 - 3d 3 / 2 lines are 
significantly weaker, although mixing among the jj configurations makes it difficult, occa- 
sionally, to unambiguously identify a line with a single-electron transition. The observable 
isonuclear trend of the 2p - 3d transition array, shifting to shorter wavelengths as the Fe ion 
charge increases, can be used to accurately diagnose the ionization range of the absorbing 
plasma. The spectral resolution of the contemporary soft X-ray grating spectrometers on 
board Chandra and XMM-Newton are sufficient for measuring the trough of the Fe 2p - 3d 
UTA rather precisely. Starting from Fe 9+ (~ 14.9 A), where the 3p subshell has one vacancy, 
and continuing to higher charge states (~ 14.0 A for Fe 15+ ) with more 3p vacancies, the 
2s - 3p lines, mostly 2s - 3p 3 / 2 , can also be seen in Fig. 3. These lines will usually not 
be detectable in the spectrum, especially since they overlap with stronger absorption lines 
that arise from 2p excitations to high-n levels (see Sec. 5). Fe 15+ has also two 2p-3s lines 
calculated at 17.395 A and 17.095 A, which are not shown in Fig. 3. 

The present atomic data were compared with results of similar, but partial, calculations 
by Raassen & Kaastra (2001) for several ions in the Fe 6+ - Fe 14+ range. Raassen & Kaastra 
(2001) used Cowan's atomic code (Cowan 1981), which employs the LS angular momentum 
coupling scheme with relativistic corrections. The mean wavelengths for individual ions 
agree to within a few mA between the two sets of results. The total oscillator strengths 
obtained with Cowan's code are systematically higher than the present results. The highest 
discrepancy of 15 % occurs for Fe 10+ , but for the most part the discrepancies are about 5 - 
10 %. These discrepancies could be explained by the different methods, but also by Raassen 
& Kaastra (2001) not including the important [3s3d + 3p 2 ] configuration mixings. 

3. SPECTRAL MODELING 

The optical depth r(z/) at frequency v is defined through the exponential attenuation 
I(y) = I (u)exp[— t(u)], where Ioiy) an d I(v) are the unabsorbed and attenuated intensities, 
respectively. 
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3.1. Individual line absorption 

The optical depth Tij{y) due to an absorption line (i — > j) can be written as: 

Tij {y) = N ion a i:j (is) (1) 

where Ni on is the ionic column density towards the source (in cm -2 ) and Oijiy) is the cross 
section for photoexcitation (in cm 2 ) from % to j. The effect of more than one line at a given 
frequency is multiplicative (a product of exponents). The photoexcitation cross section is: 

7re 2 

= TT^MM (2) 

Here, e and m e are the electron charge and mass, c is the speed of light, is the line 
oscillator strength, and <f){v) is the normalized line profile, generally described by a Voigt 
function, which is the convolution of the (Gaussian) Doppler broadening due to temperature 
and turbulence and the (Lorentzian) natural width of the line. The Doppler broadening can 
be characterized by a Doppler width Az/j> 

where \j is the line wavelength, T ion and M ion are, respectively, the temperature and mass 
of the absorbing ions, and v tU rb is the turbulent velocity. Throughout this work we use 
kTi on = 10 eV. In typical absorbing, photoionized, AGN material, turbulence broadening 
{vturb > 100 km/sec) totally dominates over temperature broadening. The natural width Fj 
(in sec -1 ) is associated with the total rate for depletion of the upper level j of the transition: 

r; = 5>*' + £^* (4) 

i'<j k'<j 

In the present work, we include both the radiative (A) and autoionization (A a ) rates for 
determining the natural width. 



3.2. Ionization Parameter Zones 

Due to the comparable ionization energies of the Fe M-shell ions and to the broadband 
nature of the ionizing flux, in physical photoionized plasmas, several charge states will always 
be present simultaneously, even if the gas is well confined in space and uniform in density and 
in temperature. Hence, the observed Fe 2p - 3d UTA will necessarily be a blend of several 
charge states. The spectral location and shape of this feature depend on the ionization 
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parameter. We adopt the definition of the ionization parameter: £ = L/n e r 2 from Tarter, 
Tucker, & Salpeter (1969), where L is the luminosity of the ionizing source, n e is the electron 
density, and r is the distance from the source. The photoabsorption spectra expected for £ 
= 1, 3, and 10 erg sec -1 cm are plotted in Fig. 4. The total Fe column density is taken to 
be 2xl0 17 cm~ 2 and the ionic fractional abundances are calculated with XSTAR (Kallman 
& Krolik 1995). The Fe 2p - 3d UTA can be seen to shift towards shorter wavelengths as £ 
increases. The dominant contributions for £ = 1, 3, and 10 erg sec -1 cm come, respectively, 
from Fe 5+ - Fe 7+ , Fe 7+ - Fe 10+ , and Fe 10+ - Fe 13+ . For £ = 10 erg sec" 1 cm, the weak 2s - 
3p lines below 15 A, which become available for Fe 9+ and higher charge states, can also 
be seen in the bottom panel of Fig. 4. The Fe 2p - 3d UTA broadens with £ due to the 
increasingly larger spacing between charge states for the high-ionization species (see Table 1 
in Sec. 4). Despite the various ions partially overlapping in the spectrum, the imprints of 
individual ions can still be noticed, even for a turbulent velocity of 500 km/sec that was 
used for Fig. 4. Each dip in the theoretical spectra in Fig. 4 can be ascribed primarily to 
one charge state. For the correspondence between (average) wavelength and charge state see 
Table 1 in Sec. 4. For single-^ systems with comparable velocities (< 700 km/sec), these 
individual ion contributions could be just marginally resolved. For a detailed analysis of the 
Fe 2p - 3d UTA in NGC 3783 and how it breaks up into single-ion contributions, see Blustin 
et al. (2001). 

Interestingly, all currently known detections of the Fe 2p - 3d UTA are more or less 
consistent with a single ionization parameter of £ = 3 erg sec -1 cm, while lower and higher 
ionization states are clearly not observed. In fact, the calculated UTA for £ = 3 erg sec -1 cm 
in Fig. 4, which falls roughly between 16 and 17 A, resembles the observed UTA in both 
NGC 3783 and IRAS 13349+2438 (see Fig. 1) fairly accurately. Particularly tantalizing is 
the absence from the soft X-ray spectra of the ions Fe 12+ - Fe 15+ , while both lower and higher 
(Fe-L) charge states show prominent absorption features. This discontinuity in ionization 
could be indicative of thermal instabilities precluding these ions from co-existing with the 
adjacent charge states. Investigating this possibility is beyond the scope of this paper, 
but preliminary calculations show that thermal instabilities could arise, depending very 
sensitively on the ionizing spectrum. The absence of charge states lower than Fe 6+ , on the 
other hand, is easier to understand, as a consequence of the typical X-ray fluxes in the AGN 
vicinity, which will instantly ionize a few electrons from the Fe atoms exposed to it. In regions 
further out from the center, where £ could be much lower, the column density through the 
outflowing and absorbing gas diminishes rapidly and, therefore, it is not observed. 
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3.3. Equivalent Width 

It is customary to assess the total amount of absorption with the equivalent width (EW), 
which is essentially the area under the (normalized) transmission curve (e.g., Figs. 2 and 3). 
The EW due to inner-shell n = 2 to 3 photoexcitation is plotted in Fig. 5 as a function of 
N ion (the so-called curve of growth) for five selected Fe M-shell ions at v turb = 500 km/sec. 
The EW behavior of these ions represents the EW trend for the other Fe M-shell ions as well. 
The EW is calculated in the 10 - 20 A range only, which is adequate because at the column 
densities considered for Fig. 5, the line wings do not digress much from that region. The EW 
can be seen to generally increase with the charge state. This is due to the gradual increase in 
the total oscillator strength with the ionic charge, which, in turn, is mostly attributed to the 
opening of more and more vacancies in the n = 3 shell. No increase in the oscillator strength 
occurs from Fe 0+ to Fe 2+ , because the ground configurations for those ions differ only by 
their 4s shell, which does not affect the n = 2 to 3 excitation. As seen in Fig. 5, saturation 
occurs at about 10 17 cm -2 , except for Fe 8+ , which has (a closed 3p-subshell and thus) only 
three 2p - 3d lines (see Fig. 3). The total oscillator strength is concentrated mostly in two 
of these lines. Consequently, they tend to saturate at lower column densities, as can be seen 
in Fig. 5. Note that at very high column densities, photoexcitation gradually plays a lesser 
role, since the absorption in photoionization edges (not included in Fig. 5) exceeds that in 
the lines by far. This effect is discussed further in Sec. 5. 



4. ABBREVIATED DATA SET 

There are, generally, several dozens of inner-shell n = 2 to 3 lines for each Fe M-shell 
ion. Typically, 20 lines from each ion account for more than 90 % of the EW in typical 
AGN conditions, with the exception of Fe 5+ , Fe 10+ , and Fe 11+ , where the 20 strongest lines 
represent a somewhat smaller fraction of that effect. As can be seen in Figs. 2 and 3, most 
of the lines for each ion are concentrated around a narrow spectral region, which can be 
very close to the corresponding regions for the adjacent charge states. Since several Fe M- 
shell ions are always present in the plasma simultaneously and since typical AGN turbulent 
velocities of a few 100 km/sec broaden the spectral lines, the entire ensemble of these lines 
pertaining to different charge states is coalesced into a broad UTA. This is independent of the 
spectral resolution of the measuring device. Consequently, the full set of the current atomic 
data, including its large number of lines for each ion, would in many cases be unnecessary. 
Alternatively, mean wavelengths and summed oscillator strengths could be adequate for 
modeling and for spectral fitting purposes. 
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Bauche-Arnoult, Bauche, & Klapisch (1985) proposed compact formulas for the mean 
wavelengths and spectral widths of transition arrays. However, since we have already cal- 
culated all of the lines level by level, we can simply use our detailed results to obtain these 
quantities. Additionally, the sum rules in Bauche- Arnoult, Bauche, & Klapisch (1985) do 
not include the effects of configuration mixing. In order to make the atomic data man- 
ageable, we have divided the lines of each charge state into a few groups. Each group is 
represented by the total oscillator strength of lines in the group E/^ and by an /-value 
averaged wavelength \ av : 

K v = (5) 

The sums are carried out over upper levels j, while the lower level i, here, is always the 
ground level. The effective natural width of a group T e ff is defined by: 

Each group also has a statistical spectral width AA defined by its /-weighted standard 
deviation: 

Assuming that the distribution of lines within a group is roughly normal (Gaussian), and 
for mere modeling purposes, the statistical width can be added to the Doppler width Aud 
as an additional "velocity" component [compare with Eq. (3)]: 




The factor 2 in front of the last term in Eq. (8) is needed to properly translate AA defined 
in Eq. (7) to the usual velocity dependence used for Av D . 

The groups of n = 2 to 3 absorption lines for Fe M-shell ions in the ground level are listed 
in Table 1. In the first three columns, the Fe ion is given together with the configuration (in 
the jj coupling scheme) and total angular momentum quantum number Jq, of the ground 
level. The n — 1 and 2 inner electronic shells not specified in the table are assumed to 
be full. Subsequently, the average wavelength A a „, the total oscillator strength S/ij, the 
effective natural width T e jf (due to autoionization and radiative decay), and the statistical 
width AA, are given for each group. The last two columns in Table 1 indicate the leading 
radiative transition in the group and the total number of (electric dipole) lines it represents. 
Both for the highly ionized and for the close-to-neutral species, the separation into groups 
nicely traces the spin-orbit array splitting discussed in Bauche- Arnoult, Bauche, & Klapisch 
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(1985), where each group can be characterized by a single jj electron transition (e.g., 2px/ 2 - 
3d 3 / 2 and 2p 3 / 2 - 3d 5 / 2 ). For intermediate charge states Fe 10+ to Fe 13+ , where the [3s3d + 3p 2 ] 
configuration mixing is strong, the actual atomic levels are often composed of components 
from more than one configuration, making the single-electron transition identification not 
absolutely strict. For those ions, the lines do not split as nicely into spin-orbit arrays, as can 
be seen in Fig. 3. 

Table 1 shows that the total oscillator strength generally increases with the ionic charge, 
while the natural width dominated by the total autoionization rates, remains more or less 
constant. Except for the lowest charge states (Fe 0+ - Fe 4+ ), the groups of strong 2p - 3d lines 
from each ion are nicely separated by more than 0.1 A, easily resolvable with Chandra and 
XMM-Newton, enabling an unambiguous measurement of the range of charge states existent 
in the absorbing plasma. Although the accuracy of the relativistic parametric potential 
method in HULLAC might slightly diminish as the charge state approaches neutral, for the 
present inner-shell transitions, these inaccuracies are expected to be much less than 0.1 A. 

As mentioned above, the statistical widths given in Table 1 facilitate the use of the 
present atomic data for spectral analysis. We compared the curves of growth obtained with 
the full atomic data set to those obtained with the abbreviated data set for several cases. 
The results agree to better than a few percent for Ni on up to 10 16 cm -2 (and then again 
for Ni on > 10 19 cm -2 ). On the other hand, since the largest statistical widths of AA ~ 
0.1 A, which occur for the strong groups of Fe 10+ - Fe 13+ , correspond roughly to velocities of 
2500 km/sec, for these ions, the abbreviated atomic data set given in Table 1 will reproduce 
the saturation of the curve of growth correctly only for similarly high velocities. For lower 
velocities, the shortened method will slightly overestimate the column density at which 
saturation takes place. In the most extreme cases we checked: Fe 10+ - Fe 13+ at 100 km/sec, 
the shortened method can overestimate the actual EW by as much as 40 % at a column 
density of 10 17 cm -2 . For lower charge states or higher velocities, the agreement is much 
better in the entire range of column densities. In short, although the grouping of lines 
can be useful for spectral fitting and modeling, the shortened method should be used with 
caution, particularly for low velocity systems (< 100 km/sec) and column densities of 10 17 - 
10 18 cm" 2 . 



5. ADDITIONAL AND RELATED ABSORPTION 

It is important to note that the present set of calculations is by no means the complete 
picture for absorption by Fe M-shell ions. The n = 2 to 3 transition array is, indeed, the 
most conspicuous in the soft X-ray spectrum, but will be regularly accompanied by weaker 
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lines emerging from excitations to higher levels. At high column densities, the photoioniza- 
tion edges could also play a major role. Additionally, even for moderately dense plasmas, 
the population of ions in excited levels, especially metastable levels, which occasionally lie 
extremely close in energy to the ground level, can be sufficient to produce photoabsorp- 
tion from excited levels. This would be due to collisional- and photo- excitation as well 
as radiative cascades towards these levels. In these cases, the absorption spectrum can be 
calculated only after the excited level populations are determined through comprehensive 
collisional-radiative models (that include photoexcitation) for each ion. 

Full models for all of the ions are beyond the scope of this paper, but we do wish to 
give an idea of the total photoabsorption effect from the ground level. For this purpose, we 
choose Fe 7+ for its relatively simple atomic structure. The absorption spectrum is calculated 
including inner-shell 2p - nd photoexcitation for n = 3 through 9, as well as the photoion- 
ization edges at the n — > oo series limit. All these are readily obtained with the HULL AC 
code. The results are presented in Fig. 6 for three ionic column densities of 3xl0 16 , 10 17 , 
and 10 18 cm -2 . The latter column density is rarely relevant to AGN sources, but nonetheless 
is used here to exemplify the high-column limit. The turbulent velocity used to calculate the 
transmission functions in Fig. 6 is fixed at 100 km/sec. At a column density of 3xl0 16 cm~ 2 , 
where the 2p - 3d lines are not yet saturated, the high-n lines can be seen to contribute to the 
absorption spectrum. In fact, the high-n lines can potentially give a sensitive indication of 
the ionic column density in the region around 10 17 cm -2 , where the 2p-3d UTA is saturated. 
The photoionization edges, however, although contributing to the EW, are spread out and 
therefore will be hard to detect. Note the unabsorbed region of the spectrum around 16 A 
between the 2p - 3d (UTA) and the 2p - 4d features. In the soft X-ray spectra of many AGN, 
but particularly in the spectrum of NGC 3783 where the Fe 2p - 3d UTA is deep (Blustin et 
al. 2001), the absence of absorption in this region creates a hump that could be mistakenly 
interpreted as emission. At very high column densities, the photoionization edges dominate 
the lines as seen in Fig. 6. Both the L-shell edges (photoionization from n = 2) in the soft 
X-ray band, but also the M-shell edges (n= 3) that occur at longer wavelengths (> 50 A), 
can have an appreciable effect on the observed spectrum between 10 and 20 A. Both these 
effects can be seen in the bottom panel of Fig. 6 in the high column density plots for N ion = 
10 18 cm' 2 . At these high column densities, even many of the high-n lines are saturated. 

In order to quantitatively compare the effect of the 2p - 3d UTA with those of the high- 
n lines as well as with those of the photoionization edges, we plot in Fig. 7 the EW as a 
function of ionic column density, separately for the 2p - 3d lines, the 2p-nd (n = 4 - 9) lines, 
and the L-shell and M-shell photoionization edges. The total EW is also shown. Similarly 
to Fig. 6, the turbulent velocity is fixed at 100 km/sec. It can be seen that up to Ni on = 
10 17 cm -2 , the 2p - 3d lines dominate the EW. At higher column densities, when the 2p - 
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3d lines saturate, both the L-shell and M-shell photoionization edges give the predominant 
absorption effect. The high-n lines give a rather small contribution to the EW, about half 
that of the 2p - 3d lines, although judging from Fig. 6, they should be detectable at ionic 
column densities of ~ 10 17 cm -2 . These numbers could vary slightly from ion to ion, but 
the general trends illustrated in Fig. 7 are representative of the EW behavior for all of the 
Fe M-shell ions. For ions, such as Fe 8+ , in which the 2p-3d lines saturate at lower column 
densities, the relative contribution of the high-n lines would be more significant. 

Note that analogous 2p - 3d UTA absorption can also be produced by elements other 
than Fe. However, in order for those to be observed, many M-shell ions pertaining to those 
elements need to show high EW in the inner-shell 2p - 3d lines, which in turn requires a 
much lower ionization parameter than that for the Fe 2p - 3d UTA. To give a general idea of 
where these lines could be expected in the spectrum, we calculated the wavelengths for the 
strongest inner-shell 2p - 3d lines of (Mg-like) Ca 8+ , Ar 6+ , S 4+ , and Si 2+ . The approximate 
wavelengths are, respectively: 31.4 A, 43.5 A, 64.0 A, and 106 A. Of these, only the Ca 2p 
- 3d UTA falls in the RGS wavelength band, while the others could be observed with the 
LETG spectrometer on board Chandra. In fact, preliminary analysis of the AGN soft X-ray 
spectra available to us show features that could, tentatively, be identified as inner-shell 2p - 
3d absorption by M-shell Ca. More details will be presented when the analysis is completed. 



6. CONCLUSIONS 

Extensive calculations for inner-shell n = 2 to 3 (mainly 2p - 3d) photoexcitation have 
been carried out for neutral iron (Fe 0+ ) and for the 15 ions Fe 1+ through Fe 15+ . The cor- 
responding absorption forms a Fe 2p - 3d UTA that has been observed in AGN spectra 
between 16 - 17 A. It is shown that despite strong blending within charge-states, the sepa- 
ration (generally > 0.1 A) between the individual ion features enables accurate diagnostics 
of the ionization range present in the absorbing medium. Ionic column densities of up to 
about 10 17 cm~ 2 can be relatively easily diagnosed. Beyond that, the Fe 2p - 3d UTA will 
be saturated and other features, such as high-n lines need to be used. Turbulent velocities 
of < 700 km/sec can still be marginally resolved. An abbreviated set of atomic parame- 
ters, which can facilitate the spectral modeling of the rather complex Fe 2p - 3d UTA, is 
provided and found in most cases to reproduce approximately the EW obtained with the 
full set of atomic data. Extra caution needs to be applied when using it for low velocities 
(< 100 km/sec) and ionic column densities of 10 17 - 10 18 cm -2 . Finally, the effects of ab- 
sorption via 2p - nd (n > 4) lines, which are natural to accompany the 2p-3d UTA, should 
be observable at ionic column densities of ~ 10 17 cm -2 . Conversely, the spectral signature 
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of the photoionization edges is harder to detect directly, although for ionic column densities 
of ~ 10 17 cm' 2 and higher, their contribution to the EW is predicted to exceed that of the 
lines. Finally, it is pointed out that 2p - 3d UTA features by lighter elements, such as Ca, 
Ar, S, and Si, could also be sought in the soft X-ray and EUV spectra of AGN. 

The authors are grateful to the rest of the RGS team for ongoing collaboration, which 
resulted in the discovery of the Fe 2p - 3d UTA feature and prompted this work. E.B. 
acknowledges useful discussions with Marcel Klapisch on the physical properties of UTAs. 
We acknowledge generous support from the National Aeronautics and Space Administration. 
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Table 1 — Continued 
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a total angular momentum quantum number of the ground level. 
b denned in Eq. (5). 
c denned in Eq. (6). 
d denned in Eq. (7). 
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Fig. 1. — The RGS spectra of the Seyfert 1 NGC 3783 and the luminous Quasar IRAS 
13349+2438 featuring a broad trough between 16 - 17 A, due to numerous inner-shell n— 2 
to 3 (mainly 2p - 3d) absorption lines of Fe M-shell ions blended into a UTA. 
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Fig. 2. — Spectra of inner-shell n = 2 to 3 absorption lines for Fe 0+ through Fe 7+ . A 
temperature of kT ion = 10 eV and no turbulent velocity are assumed for these plots, in order 
to minimize line blending and emphasize the complex atomic structure and rich line spectra. 
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Fig. 3. — Spectra of inner-shell n = 2 to 3 absorption lines for Fe 8+ through Fe 15+ . A 
temperature of kT ion = 10 eV and no turbulent velocity are assumed for these plots, in order 
to minimize line blending and emphasize the complex atomic structure and rich line spectra. 
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Fig. 4. — Multi-ion spectra featuring numerous inner-shell n = 2 to 3 (mainly 2p - 3d) 
absorption lines of Fe M-shell ions blended into a UTA. Calculations are carried out for 
single ionization-parameter zones (£ = L/n e r 2 ). N Fe = 2xl0 17 cm~ 2 and v tur b = 500 km/sec 
are assumed for the three plots. 
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Fig. 5. — Curves of growth for five selected Fe M-shell ions due to the inner-shell n = 2 to 
3 (mainly 2p - 3d) absorption lines. EW is calculated in the wavelength region between 10 
and 20 A. v tU rb is taken to be 500 km/sec. 
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Fig. 6. — Spectra of Fe 7+ at three different ionic column densities. Calculations include all 
of the inner-shell 2p - nd (n = 3 to 9) absorption lines, as well as the L-shell and M-shell 
photoionization edges. The turbulent velocity is fixed at 100 km/sec. 
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Fig. 7. — Curves of growth for Fe 7+ . Separate curves are given for the EW contributions of 
inner-shell 2p - 3d absorption lines, inner-shell 2p - nd (n = 4 - 9) absorption lines, and for 
the L-shell and M-shell photoionization edges. EW is calculated in the wavelength region 
between 10 and 20 A. v tur b is taken to be 100 km/sec. 



